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Studies of electronic charge transport through semiconductor double quantum dots rely on a 
conventional 'hole' model of transport in the three-electron regime. We show that experimental 
measurements of charge transport through a Si double quantum dot in this regime cannot be fully 
explained using the conventional picture. Using a Hartree-Fock (HF) formalism and relevant HF 
energy parameters extracted from transport data in the multiple-electron regime, we identify a novel 
spin-flip cotunneling process that lifts a singlet blockade. 
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In quantum computing, semiconductor quantum dots 
have long been considered as good candidates for 
qubits [TH3]. A promising architecture for such qubits 
is the double quantum dot [SHE]. Understanding spin- 
dependent transport [6UTU] is important for using the 
spin degree of freedom in a double dot qubit. Here, 
we show that transport data taken in the three-electron 
regime of a double dot in a Si/SiGe heterostructure have 
features that are qualitatively inconsistent with the con- 
ventional model of 'hole' transport [11 , because this 
model does not account for transport through excited 
states. Using the Hartree-Fock (HF) formalism with 
singly excited configurations [12] . together with relevant 
HF parameters extracted from the transport data (see 
[13]), we demonstrate that the striking features in the 
data arise from a novel spin-flip cotunneling process in 
which the multi-electron nature of the system enters fun- 
damentally. 

Several experiments have probed charge transport 
through double quantum dots in the few-electron regime 
and investigated effects such as energy-dependent tun- 
neling and spin-dependent transport pHTU]. Transport 
in the three-electron regime is well-described in terms of 
holes when all the intra-dot relaxation rates are much 
faster than the interdot tunnel rate, so that the domi- 
nant transport channels are through the lowest energy 
states of each dot, as is typically the case in GaAs de- 
vices [1 EE] . 

Our theoretical work is based on data [8 , in which a 
lateral double quantum dot was formed by electrostatic 
gating of a Si/SiGe heterostructure, as shown in the inset 
of Fig. [TJa). Fig. [TJa) shows source-drain current ver- 
sus controlling gate voltages at a fixed source-drain bias 
voltage. Transport through the two dots is energetically 
favorable within triangular regions whose size is deter- 
mined by the source-drain bias. Lines of high current in 
these bias triangles are associated with fast tunneling be- 
tween the dots and between the dots and the leads [TT] . 
From the orientation of the line a/3 in Fig. [ljb) , we de- 
duce that it is associated with the resonance of an energy 
level in the 'left' dot with the chemical potential in the 
left lead. Quantitative fits allow the edges of the triangles 



to be determined and are reported in detail in Ref. [13] 
Fig. [ljb) is a schematic diagram of the bias triangles, 
with energy axes shown in the inset. The lower features 
arise from transport when the double dot contains either 
one or two-electrons ('two-electron' regime), while the 
upper features reflect transport when the dot contains ei- 
ther two or three-electrons ('three-electron' regime, also 
conventionally termed 'hole' regime). 

There are two regions of current flow in each transport 
regime, shown in Fig. [TJa). Each of these regions of cur- 
rent is contained in a triangle, shown in either blue or red 
in Fig.JTJb). The presence of current in the blue triangle 
implies that there is significant transport through excited 
states of the dots [8j [13], something that has recently 
been observed in transport through a single phosphorous 
donor in silicon as well [T4] . 

Because the effective electronic mass in Si is much 
larger than in GaAs, transport is energetically favor- 
able within each bias triangle, but the triangle is not en- 
tirely filled because the electron tunneling rate is strongly 
energy-dependent [6H5]. The two parallel lines of high 
current that are observed along the left edges of the sin- 
glet and triplet triangles in the two-electron regime (lower 
feature) indicate that energy-dependent tunneling across 
the left barrier is the bottleneck in the total tunneling 
rate [T3] . 

In the three-electron regime, the conventional picture 
that describes the conduction in terms of holes predicts 
that there should be two parallel lines of high current 
(Fig. [TJc)). This is because, in the two-electron regime, 
electron occupancy cycles through the states (1,0) — >• 
(2,0) — >• (1,1), and the three-electron regime is modeled 
conventionally [TT] as hole transport in the opposite di- 
rection: (1,1) — ^ (0,2) — )> (0,1), where the numbers repre- 
sent electron or hole occupancy in the left and right dots. 
Due to particle-hole symmetry, the hole picture predicts 
parallel lines of high current similar to the data in the 
two-electron regime (Fig. [TJc)). 

The transport data in Fig. [TJa) is inconsistent with a 
picture in terms of holes, as it shows two lines of high 
current in the three-electron regime (upper feature) that 
are clearly not parallel. In this regime, there is a line of 
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FIG. 1: (a) Transport current /sd in a Si/SiGe double quantum dot (color scale) as a function of controlling gate voltages, Vb 
(V) and Vcs (V), reported in Ref.[8| Inset shows an SEM image of the gates with a numerically simulated double dot overlayed. 
White letters (T, Bl, G, Br and CS) label gates and red letters (S, D) label the source and drain. For this data, electrons flow 
from left to right, (b) A cartoon of the bias triangles and lines of high current. Inset shows the energy axes of the dots. The 
lower (upper) features are the electron (hole) triangles. Red dashed lines represent current through the singlet (singlet-like) 
channel of the electron (hole) bias triangle and blue dot-dashed lines the triplet (triplet-like) triangle. A and B are resonant 
peaks of the singlet and triplet electron triangles. C is the resonant peak of the singlet-like hole triangle. A and C are the 
triple points at the boundary of the (1,0), (1,1), (2,0) and (1,1), (2,0), (2,1) charge occupations. F lies along the line extending 
from the tail; it is a representative point where cotunneling is dominant. Est is the (2,0) singlet-triplet energy splitting. Data 
are obtained at a reverse-bias source-drain voltage, Vsd = —0.274 mV, first published in Ref. 8 as —0.3 mV. Ref. 1131 details the 
quantitative fits to identify the triangles, (c) The prediction using the conventional hole picture in the three-electron regime is 
shown as two parallel lines (black), which is inconsistent with the tail observed in the data. 



high current at the left edge of the bias triangle (line £77) 
for ground state transport, which is expected since the 
left barrier is observed to be the bottleneck in the two- 
electron regime. However, in the bias triangle for excited 
state transport, there is a 'tail' parallel to the right edge 
but away from it, which the hole picture completely fails 
to describe. 

To understand the problem theoretically, we formulate 
it in terms of chemical potentials and use the Hartree- 
Fock (HF) approximation with singly excited configura- 
tions to determine the spin eigenfunctions and energy 
levels of each of the double dot states involved in the 
three-electron regime. The relevant parameters in the 
HF formulation are extracted from the transport data as 
detailed in Ref. [13l From the energy levels and possible 
transitions between states, we calculate the electrochem- 
ical potentials for charging or discharging a dot by one 
electron [15 j. The four relevant electrochemical poten- 
tials for the dots are shown for the two-electron case in 
Fig. [2ja,b). For the three-electron case, the full set of 
ten electrochemical potentials, shown in Fig. [2|c,d), is 
clearly greater than the four electrochemical potentials 
for transport modeled on two holes. The many-electron 
nature of the problem thus enters our analysis of trans- 
port naturally. 

Without going into the details of the HF calculations, 
we can gain some insight into the possible (2,1) states us- 
ing qualitative arguments. The pure singlet and triplet 
states, 5(2, 0) and T(2, 0), are no longer orthogonal when 



we include a weak coupling to a third electron in the right 
dot. The perturbation leads to a 'singlet-like' ground 
state 5* (2, 1), whose spin configuration in the left dot is 
mainly 5(2,0) with a small admixture of T(2,0). The 
5* (2,1) state has spin S z = ±1/2 and is doubly de- 
generate. The perturbation also leads to 'triplet-like' 
states T*(2, 1), for which spin addition gives S z = ±1/2 
or ±3/2. The S z = ±1/2 states contain mainly triplet 
T(2, 0) with a small admixture of 5(2, 0). The S z = ±3/2 
states have spins that are either all up or all down; they 
are doubly degenerate without any admixture of singlet 
states. The triplet degeneracies are lifted due to the fact 
that exchange energies are different for different three- 
electron spin configurations. The energy splittings arise 
from inter-dot interactions, which are much smaller than 
intra-dot interactions. Thus, the splittings within the 
triplet-like manifold are much finer than the splitting be- 
tween the singlet- and triplet-like manifolds. These ar- 
guments are borne out by our calculations [T5] . 

From the energy levels calculated with the HF Hamil- 
tonian, we can explain how the electrochemical po- 
tentials, shown in Fig. |2|c), are obtained. In the 
three-electron regime, electron occupancy cycles through 
(1, 1) — ^ (2, 1) — >> (2,0). The first transition corresponds 
to charging of the left dot from a (1,1) to a (2,1) state. 
For clarity, we do not distinguish between the two closely 
spaced (1,1) energies, nor do we distinguish between the 
three closely spaced T*(2, 1) energies. We therefore ob- 
tain two distinct electrochemical potentials, /jL c ,t* an d 
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FIG. 2: Diagrams for transport through excited states and 
the process of 'spin-flip cotunneling.' (a) Triplet channel 
transport in the two-electron regime. At B, the resonant peak 
of the triplet channel, transport is allowed through the triplet 
levels, /xt(2,o) and /4t(i,i)- (b) When the singlet level /is(2,o) 
is loaded, transport is energetically uphill and blockaded. Co- 
tunneling of a left dot electron out to the right lead lifts the 
blockade to resume triplet channel transport, (c) In the three- 
electron regime, transport occurs in the following cycle: step 
1, (1, 1) (2, 1); step 2, (2, 1) (2, 0); step 3, (2, 0) (1, 1). 
At D, these occur through the triplet- like channel (blue), (d) 
Because of the loading of the singlet-like /i c ,s* level (red) at 
D, the system ends up in the 5(2, 0) state, whereby transport 
is blockaded, (e) When an electron cotunnels from the left 
lead into the right dot to form a triplet-like state, it puts the 
left dot into an admixture of singlet and triplet, (f) The right 
dot discharges from fid,T* , leaving a triplet state in the left 
dot, thus causing a spin- flip and resuming transport. 



Hc,s* •> shown in Fig.[2^c), which are the energies needed to 
charge the left dot from a (1,1) state to the T*(2, 1) and 
5* (2, 1) states respectively. The second transition repre- 
sents the discharge of an electron from a (2,1) to a (2,0) 
state. Electrochemical potentials, Hd,T* and ^d,s* •> drawn 
on the right dot, represent the discharge of the right dot 



from T*(2,l) to T(2,0) and 5* (2,1) to 5(2,0) respec- 
tively. These are the continuous lines (blue and red) on 
the right dot in Fig. [5|c). Due to singlet-triplet mix- 
ing in the left dot, two other transitions of much smaller 
likelihood are possible. They are the 5* (2,1) to T(2,0) 
and T*(2, 1) to 5(2, 0) transitions, represented by the red 
dotted and blue dashed levels respectively, in the same 
figures. The last step in the cycle is the inter-dot transi- 
tion, (2,0) — >• (1, 1). The chemical potentials in this step 
are identical to the two-electron case (Fig. [2|a,b)) and 
are labeled as /is ? T(2,o) an d Ms,t(i,i)- 

We can now explain the tail in the transport data, 
which, as described above, is a prominent feature that 
is qualitatively inconsistent with a description in terms 
of holes. At point D in Fig. [2|c), transport is allowed 
through the blue triplet-like levels. However, it is also 
possible to load the red singlet-like ji C: s* level. In this 
case, as the right dot discharges, the system is likely to 
end up in the 5(2,0) state, where transport is energeti- 
cally uphill and therefore blockaded (Fig. [2jd)). We call 
this a 'singlet blockade.' The lifting of the singlet block- 
ade along the tail is shown in sequence in Figs. [2je) and 
(f). Starting from 5(2, 0), the double dot forms a triplet- 
like T*(2, 1) state when an electron from the left lead 
cotunnels into the right dot, as shown inFig.^e). The 
charging of the right dot in this transition requires the 
same energy as its reverse discharging process (T*(2, 1) to 
5(2, 0)), represented by the blue dashed line in Fig.[2|c). 
It is labeled by /i C 2,T* in Figs. [2|e) and (f). Because 
the triplet-like state contains an admixture of singlet and 
triplet states in the left dot, when the right dot discharges 
from the jid,T* level, the left dot ends up in the triplet 
(2,0) state, thus causing a spin-flip. With the singlet 
blockade lifted, the system then completes the cycle into 
the triplet (1,1) state and transport resumes as shown in 
Fig. [5|c). We term this process 'spin- flip cotunneling.' 

The tail in the transport data in Fig. [TJa) is bright 
along its entire length because the chemical potentials 
for the right dot and the left lead are the same. Point D 
is the brightest point along the tail because of the fast 
inter-dot tunneling when /iT(2,o) I s aligned with /iT(i,i)- 

The spacing of the tail away from the edge of the tri- 
angle is consistent with the energy difference between the 
Mc2,t* and /id,T* levels on the right dot (Fig.[2|e)) being 
equal to Est, the (2,0) singlet-triplet energy splitting. 
To understand how this is consistent with the transport 
data, we start from point C in Fig. [TJb) and note that 
when both dot energies fall by Est, the blue, dashed 
Mc2,t* level of the right dot lines up with the Fermi level 
of the left lead. This measure of Est is also consistent 
with other measures of ST splitting [13] . 

The significant role cotunneling plays in the triplet 
and triplet-like transport channels of the two and three- 
electron regime is interesting. In both cases, cotunneling 
by itself does not contribute significantly to the current, 
but plays the role of allowing transport to resume by 
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lifting the singlet blockade. 

Current will flow through the triplet-like channel when 
the loading rate is comparable to the unloading rate in 
the singlet-like channel [15] . In the Supplementary Infor- 
mation, we estimate these rates and find that they are 
indeed the same order of magnitude. The blockade is 
therefore lifted about as quickly as it is encountered. In 
this way, spin-flip cotunneling enables transport through 
the triplet-like channel. The resulting current is that of 
the unblockaded, triplet-like channel, reduced by a factor 
of - 2 [15]. 

Interestingly, transport via the triplet channel was not 
observed in the experiments reported in Ref. [9] In that 
study, a conventional hole model was sufficient to de- 
scribe transport in the three-electron regime, as consis- 
tent with the fact that transport occurred through the 
ground states in the two-electron regime. 

It is also interesting to compare the intra-dot spin-flip 
times with inter-dot tunneling times for GaAs and Si. In 
GaAs devices, spin-flip times range from —200 /is for a 
two-electron dot [16], to —0.85 ms (at 8 T [17]) and > 1 s 
(at 1 T and 120 mK [18 ) for single electron dots. Re- 
cent experiments report spin-flip times in single electron 
dots in Si [T9H22] ranging from 40 ms (at 2 T [19]) to 
6 s (at 1 T [20]), at low temperatures. The tunnel cou- 
pling for the same Si double dot studied here was found 
to be 10 ns (25 ns) in the elastic (inelastic) tunneling 
regime [13]. However, tunnel couplings for electrostati- 
cally gated semiconductor double dots are tunable and 
can be both larger or smaller than spin-flip times. 

We note that a pulsed gate experiment [7] exhibit- 
ing phenomena arising from singlet-triplet mixing, as de- 
scribed above, is presented in [T5] . 

In summary, we have shown that the conventional hole 
model of transport in the three-electron regime fails qual- 
itatively because of the importance of excited state trans- 
port. The Hartree-Fock formalism, with relevant param- 
eters fitted to transport data, leads to the description of 
a model which explains all of the features of the transport 
data, including a novel process of spin-flip cotunneling. 

We thank Nakul Shaji for useful discussions and 
gratefully acknowledge funding from ARO and LPS 
(W911NF-08-1-0482) and NSF (DMR-0805045). 

Supplementary Material: Unconventional Transport 
in the 'Hole' Regime of a Si Double Quantum Dot 

The following Supplementary Information presents our 
detailed calculations of the Hartree-Fock approximation 
and electrochemical potentials relevant for reverse bias 
transport in the three electron regime. We estimate 
the spin-flip cotunneling occurrence rate, quantify the 
effect spin-flip cotunneling has on transport current in 
the triplet-like triangle, and show that it is consistent 
with the data. We propose an experiment to verify the 



singlet-triplet mixing mechanism and also present a de- 
tailed analysis showing that at forward bias, transport 
data has a feature associated with the process of spin- 
flip cotunneling, as discussed in the main paper. 

Hartree-Fock Approximation with Singly Excited 
Configurat ions 

The Hartree-Fock (HF) Hamiltonian used in our cal- 
culations is given by [12] 

H = ^\ h \^) a \ a j + \ Y {ij\kl)a\a]aia k , (SI) 

where the sums run over the spin orbitals chosen to 
form the basis of the problem. Operators a\ and cre- 
ate and annihilate an electron in the z-th spin-orbital 
respectively. The first term contains the one electron 
operator, h = J]— 1/2 V^, + V(x a — Xdot), the sum 
of the kinetic and potential energies of each a-th non- 
interacting electron in a quantum dot centered at x*dot- 
We note that the parabolic well approximation is a rea- 
sonable form of the potential for our geometry. The 
second term contains the two electron Coulomb opera- 
tor, (ij\kl) = J dx a dx /3 x*(x a )x*(x /3 )r^x/c(x a )x^(x / 3), 
where Xii^a) denotes the z-th spin-orbital wave- function 
of the a-th electron, following the notation of Ref. [12j 

We define the quantities sij = the overlap be- 

tween the i-th and j-th spin orbitals; e$ = (i\h\i), the 
single particle energy; Jij = and = the 

Coulomb and exchange interactions respectively, = 
(i\h\j) and Y\ 3 k = (ij\kl) govern coherent tunneling be- 
tween dots and are first and second order in the overlap 
integral, respectively. 

Singly-excited configurations that form the basis for 
the (2,1) Hamiltonian matrix, are constructed from dot- 
centered single particle spatial orbitals {1,2} and {3,4}, 
centered on the left and right dots respectively. The two 
sets of spatial orbitals are not orthogonal to each other, 
due to the overlap across the finite barrier between the 
dots. Within each set, however, the spatial orbitals are 
orthonormal. The corresponding spin orbitals are de- 
noted by {1, 1, 2, 2} and {3, 3, 4, 4} where the overbar (or 
lack of) denote spin down (up). Table [i] shows the basis 
states used to construct the 10 x 10 Hamiltonian matrix. 

The (2,1) Hamiltonian consists of four block diago- 
nals: two lxl blocks corresponding to a total spin 
of S z = ±3/2, and two 4x4 blocks corresponding to 
S z = ±1/2. In the S z = ±1/2 blocks, off-diagonal ma- 
trix elements are second order in the overlap between 
spin orbitals centered on different dots and represent the 
amount of mixing of the singlet and triplet states in the 
left dot. As such, we term the ground state, 'singlet- 
like', and denote it as S* (2, 1). It contains mainly singlet 
5(2,0) with some admixture of triplet T(2,0) in the left 
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TABLE I: Basis states used for the (2,1 
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dot, coupled to the third electron in the right dot. The 
first excited state is the pure S z — ±3/2 state with all 
spins either up or down, which does not have any singlet- 



triplet mixing. The next two higher states, which are 
closely spaced with each other as well as with the first 
excited state, contain a small admixture of singlet with 
mainly triplet states in the left dot, coupled to the right 
dot. For brevity, we term the three excited states triplet- 
like T*(2, 1) states, and do not distinguish between them 
henceforth, as the energy separation between the triplet- 
like states is of second order in the overlap integral and is 
small compared to the spacing between the ground and 
first excited state. The next higher state is well separated 
from these and does not enter into our analysis. 

The S z = 1/2 block of the Hartree-Fock Hamiltonian 
for the (2,1) state in the basis |S (2, 0)) <g> |3), |T (2,0))<g> 
|3), |Si(2,0)> ® |3),|T + (2,0))®|3) is 



/ 2ei + e 3 + Jn + 2 Jis 
— K13 + 2si 3 ti 3 
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FIG. 3: Electrochemical potentials \i for different states rele- 
vant to transport in the three electron regime. Here, the label 
Est is the (2,0) singlet-triplet splitting. A is the spacing be- 
tween the chemical potentials of the two and three electron 
regimes, as defined in Eq. |Sll| and is equal to the energy spac- 
ing between points A and C in the transport data in Fig. 1(b) 
of the main text. Both quantities are extracted from the 
transport data as detailed in Ref. 1131 giving Est ~ 174 fieV 
for (2,0) singlet-triplet splitting and A « 500 /xeV. 

The off-diagonal elements are second order in the over- 
lap integral and therefore small, so the problem can be 
treated perturbatively. The energies, to first order ap- 
proximation, are given in Table [TTJ The splitting between 
the singlet-like and triplet-like states is approximately 
equal to the singlet-triplet splitting of the (2,0) state, 



Est(2,q)- We apply the HF Hamiltonian (Eq. SI) to cal- 
culate the singlet and triplet energies of the two electron 
states, which are necessary for the calculation of elec- 
trochemical potentials. Similar analyses using molecular 
and dot-centered orbitals have been done elsewhere [23] , 
so we do not reproduce the calculations here. We refer 
the reader to Section III. D of Ref. [23] for calculations 
using dot-centered orbitals. We caution the reader that, 
in this previous study, the double dot geometry and ma- 
terial are different, so their numerical results will not be 
applicable for our case. However, the expressions for en- 
ergy levels in terms of dot-centered single-particle, direct 
and exchange Coulomb energies in Eq. (18-21) of Ref. [23] 
are relevant for our two electron double dot case. 

Once the expressions for the energies are known, it is 
straightforward to calculate the relevant expressions for 
electrochemical potentials. For clarity, we neglect the 
singlet-triplet splitting in the (1,1) states and the energy 
spacing between the three triplet-like (2,1) states and ob- 
tain two distinct electrochemical potentials. The electro- 
chemical potentials for charging the left dot from a (1,1) 
state to the S*(2, 1) and T*(2, 1) states are, respectively, 

/J>c,s* = ei + Jn + J13 - E S t (1,1) + E LdoU (S3) 

Mc,T* ~ Mc,^* + E ST (2,0)- (S4) 
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TABLE II: Energies and characteristics of the (2,1) states. 



Energy (ascending order) 




State 


Spin S z 


2ei ± e 3 ± Jn ± 2Ji 3 
-K13 ± 2si 3 £i 3 




Singlet-like, 
iviixea 


±1/2 


ei + e2 + e3 + J12 + J13 
+J23 - (#12 + #13 + #23) 




Triplet- like, 
Pure 


±3/2 


ei + e2 + e3 + J12 + J13 ± J23 — K12 

-(Kl3 + i^ 23 )/2 - (Sistis + S23t 23 ) 


Triplet- like, 
Mixed 


±1/2 


ei + e 2 + e 3 + J12 + J13 + J23 


-i^!2 


Triplet-like, 
Mixed 


±1/2 



ERdot (ELdot) is the electrostatic energy of the right (left) 
quantum dot. The electrochemical potentials for dis- 
charging the right dot from the triplet-like T*(2,l) to 
the triplet T(2, 0) state and from the singlet-like 5* (2, 1) 
to singlet £(2,0) state are 



Md,T* ~ l^d,s* = 63 ± 2Ji 3 - K 13 ± E Rdot . 



(S5) 



There are two more electrochemical potentials involved 
in the discharge of the right dot from the (2,1) states, 
due to the mixing of the singlet and triplet (2,0) states. 
These are associated with the discharge from T*(2, 1) to 
5(2, 0) and 5* (2, 1) to T(2, 0), and are shown by the blue 
dashed and red dotted lines on the right dot in Fig. [3] 
They are separated from jid,T* and jid,s* by the (2,0) 
singlet-triplet splitting Est- 

The last set of electrochemical potentials needed are 
those for the two electron regime, which are given by 

Ms(2,o) = e i + J 11 + E LdoU (S6) 

MT(2,0) = /4S(2,0) + E S T(2,0), (S7) 

Mt(i,i) = 63 ± J13 - K13 ± E RdoU (S8) 

= MT(1,1) - ^5T(1,1)- (S9) 

The chemical potentials for the two and three electron 
regimes satisfy the relation 

Mc,S* _ MS(2,0) ~ A^S* _ MS(1,1)> (S10) 

and we therefore define 



- MS(1,1) 



(Sll) 



Here, we point out that it is the relative positions of 
the electrochemical potentials that are important in or- 
der to build the transport model shown in Fig. [3j The 
important energy quantities are therefore, the (2,0) and 
(1,1) singlet-triplet energy splittings E S t(2,o) ~ 174 fieV 
and E$t(i,i) ~ 4 /ieV, and the spacing A « 500 /xeV 
between the two and three electron triangles. These en- 
ergies are extracted from the transport data as detailed 
in Ref. US 



Effect of Spin-Flip Cotunneling on Transport 
Current 

In this section, we show that current can flow through 
the triplet channel in the 'hole' transport regime when 
certain conditions are met for the singlet channel. Specif- 
ically, the singlet unloading rate should be equal to or 
larger than the singlet loading rate. We then show that 
these conditions are met in our data. 



Conditions on the singlet tunnel rates 

We first obtain the conditions on the singlet transport 
rates by developing rate equations, in analogy with those 
found in Ref. E 

We consider the contribution of loading and unloading 
rates of each channel to the total transport rate through 
the channel. Starting from the (1,1) state, the double dot 
is loaded into one of the (2,1) states from the left lead. 
The mean loading rate is given by 



Tload = l/^load 

= load ± Tt load, 



(S17) 



where T S /t load are the loading rates into the singlet-like 
and triplet-like states from the left lead respectively. 
The total time required for a complete transport cycle 

is 



T = Ti oad ±p T /r T 

unload + Ps/T s unload 1 



(Sis) 



where T$/ T unload are the unloading rates of the singlet- 
like and triplet-like channels respectively, and 



Ps 
Pt 



S load 



load ± Tt load ' 
Tt load 



■ S load 



Yj 1 



load 



(S19) 
(S20) 



are the probabilities of loading the singlet-like and 
triplet-like channels. 

In the triplet-like transport regime, the singlet-like 
channel is strongly suppressed by energy dependent tun- 
neling. If we assume that the loading and unloading rates 
for the triplet-like channel are much faster than the load- 
ing and unloading rates for the singlet-like channel, and 
if we also assume that the unloading rate of the triplet- 
like channel is much faster than its loading rate due to 
asymmetric tunnel barriers (see Ref. [8| for a discussion) , 
then Eq. S(18) reduces to 



J- S load 
r S unload 

load 
r S unload 



1 



1 



Tt load 
1 

Tt load 



T unload 



(S21) 
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In the limit of only one fast triplet-like channel, the 
total tunneling time is T « 1/IY load- This implies that if 
Ts load <^ unload, then the current through the double 
dot is essentially a triplet-like channel current. Therefore, 
it is important to determine the ratio of the loading to 
unloading rates of the singlet-like channel. 



Loading rates 

We now estimate the loading rates for the singlet- 
like and triplet-like channels at point D on the 'tail' of 
the transport data (Fig. 1, main text). Starting from 
the (1,1) charge configuration, the loading rates into the 
(2,1) states are strongly energy-dependent at the left bar- 
rier 13 and are given by 



-E S / T /E L 



(S22) 
the energy 



1 S/T, load — 1 0<5 

where the amplitude r « 1.5 x 10 8 s~ 
dependent coefficient El ~ 40 /xeV, Et 
174 /xeV [13] and Es ~ 2E$t(2,o) a ^ point D. 

The loading rates of the triplet-like and singlet-like 
states at point D are therefore 



E 



ST (2,0) 



IV, load « 1.9 XlO 6 S-\ 

T s , w « 2.5 xlO 4 s- 1 . 
Unloading rates 



(S23) 
(S24) 



We now estimate the unloading rates for the singlet- 
like and triplet-like channels at point D. Unloading of the 
triplet-like (2,1) state occurs when the electron in the 
right dot tunnels out across the right barrier, followed 
by interdot tunneling from the (2,0) to the (1,1) state. 
The tunneling rate across the right barrier T# ^ 9.2 x 
10 9 s _1 and the coherent interdot tunneling rate I\ « 
7.7 x 10 8 s _1 [13 j. The unloading rate for the triplet-like 
channel is therefore given by 



T, unload 



(r- 1 + 17 1 )- 1 « 7.1 x 10 8 



(S25) 



The unloading rate of the singlet-like state is domi- 
nated by spin-flip cotunneling, which involves two low- 
probability processes: the spin-flip and the cotunneling. 
The spin-flip can be understood as the small probability 
p m that the T* state contains the (2,0) singlet. We then 
have 



S, unload 



(S26) 



where Tc is the total cotunneling rate. 

We can estimate p m from the (2,1) Hamiltonian of 
Eq. (S2) by considering the subspace spanned by the 
states 150(2,0)) <g> |3) and |T (2,0)> <g> |3). The Hamil- 
tonian within this subspace is given by the top-leftmost 



2x2 matrix block in the right-hand side of Eq. (S2). The 
off-diagonal element in this subspace is (r^ + ^13^23 + 
^23^i3)/v // 2 « The singlet-triplet mixing proba- 

bility is therefore 



Pn 



pl3 
1 23 



2E 



ST (2,0) 



(S27) 



We do not have a direct experimental estimate for the 
quantity T^. However, we can estimate it using the 
parabolic well approximation, where the eigenstates of 
each dot are those of the 2-dimensional harmonic oscil- 
lator [24]. For a dot size of approximately 100 nm and 
an interdot distance of approximately 200 nm [8] , we ob- 
tain r3>3 ~ 0.15 meV when there is maximum overlap be- 
tween the first excited eigenstate ('2') of the left dot with 
the ground eigenstate ('3') of the right dot. This corre- 
sponds to the axis of the wavefunction of the first excited 
eigenstate being aligned with the axis of the double dot. 
A more realistic estimate where the angle between the 
wavefunction axis and the axis of the double dot is 45 
degrees, gives 



F 13 
1 23 



0.073 meV. 



(S28) 



Given that E S t (2,0) 7 
mixing probability is 



174 peV [13 , the singlet-triplet 



0.044. 



(S29) 



The cotunneling rate Tc from the left lead into the 
right dot and out to the right lead, can be estimated 
from point F of the transport data (Fig. S2 and Fig. 1(b) 
of the main text). This cotunneling has contributions 
from both the singlet and triplet channels. The singlet- 
like channels, 5, T(2,0) — » 5* (2,1), are suppressed by 
energy dependent tunneling, as is one of the triplet- 
like channels, T(2,0) T*(2, 1). The other triplet- 
like channel, 5(2,0) -> T*(2,l) is suppressed by the 
slow spin-flip process. In principle, we cannot separate 
these contributions at point F. However, we can consider 
a worst case scenario where the cotunneling current at 
point F is dominated by the channels without spin-flip, 
5(2, 0) 5* (2, 1) and T(2, 0) T*(2, 1), with very lit- 
tle energy dependent tunneling. In this case, the co- 
tunneling rate through the triplet-like spin-flip channel, 
5(2,0) T*(2,l), is given by Eq. (S26), where eT c is 
the total current measured at point F. In this way, we 
obtain 

Tc^S.OxlO^- 1 , (S30) 
T s , unload^ 2.2 xloS" 1 . (S31) 

Conclusions 

It is clear from the preceding analysis that the loading 
and unloading rates satisfy the conditions for Eq. (S21) 
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Mc2,T* 














MT(2,0) 


MT(1,1) 


MS(1,1) 






MS(2,0) 





Mi? 



FIG. 4: Electrochemical potentials at point F of the transport 
data shown in Fig. 1 of the main text. The dominant current 
contribution is from cotunneling from the left lead into the 
right dot and out to the right lead. 



to be valid, i.e. T T unload > IV load > load / unload . 
Since the ratio Ts lo&d/^s unload ~ 1? the current along 
the tail is of the same order of magnitude as the current 
in the singlet-like triangle away from the resonant left 
edge. 

Below the tail, cotunneling is energetically unfavor- 
able, while above the tail, cotunneling is suppressed 
through energy dependent tunneling. Ts unload is there- 
fore much smaller on both sides of the tail. Loading of 
the singlet-like channel, which is strongly energy depen- 
dent, increases below the tail, further suppressing cur- 
rent. However, suppression of the singlet-like channel 
loading rate above the tail could mitigate the decrease 
in the unloading rate, so that the ratio Ts \oad/^s unload 
remains small, allowing some current flow. That there is 
small current above the tail is observed in transport data 
at larger biases [13], where the loading of the singlet- 
like channel can be more strongly suppressed as it moves 
lower in the bias window. 

Therefore, although spin-flip cotunneling involves two 
mechanisms of lower probability, namely cotunneling and 
triplet-singlet mixing, the dampening effect on the tunnel 
rate is mitigated by the small loading probability into the 
singlet-like channel. The effect of spin-flip cotunneling is 
to open a thin window in the triplet-like bias triangle to 
allow triplet-like channel current. 



Proposed Experiment to Verify Singlet-Triplet 
Mixing Mechanism 

In this section, we propose an experiment to verify the 
singlet-triplet mixing mechanism described in the main 
text. We note that in the spin- flip process described in 
the text, the loading of the third electron into the right 
dot is through cotunneling from the left lead. Since it is 
the singlet-triplet mixing mechanism that we are inter- 
ested in, this loading can also be from the tunneling of 




NC r \ * 

x v « (2 11 
X % 1 ) 


7 b 






MP J 
/ 


(2,0) 


0.0) 





A*T(2,0) /id ' T * 
— Vd,S* 

/AS(2,0) 

MT(1,1) 

^(1,1) 



FIG. 5: Proposed experimental verification of the singlet- 
triplet spin-flip process. The boxed panel shows a schematic 
of the honeycomb charge stability diagram [11] , showing the 
regions in which the ground (singlet and singlet-like) states of 
the various charge configurations are stable. Points a, b and 
c are the positions to which the double dot should be pulsed. 
The dotted line indicates the boundary between the (1,1) and 
triplet-like (2,1) regions, i.e. within the shaded region, a (1,1) 
state is more stable than the triplet-like (2,1) state. Panels a, 
b and c show the electrochemical potentials, with the circles at 
the bottom of each diagram indicating the number of electrons 
at the start of each pulse. Black arrows indicate the tunneling 
of an electron, which can be detected via charge sensing using 
an adjacent quantum point contact. Detection of a tunneling 
out process at 'c' would indicate a singlet-triplet spin-flip. 



an electron from the right lead, as described below. 

Using pulsed gate techniques such as those described in 
Ref. [7J the double dot could be first loaded into the (2,0) 
singlet state (a), then pulsed such that the fi C 2,T* level 
of the right dot is below the Fermi level of the right lead 
(b), allowing for the loading of an electron into the right 
dot from the right lead. Loading of the /jl^s* l eve l will 
not cause singlet-triplet mixing but this is suppressed be- 
cause of energy dependent tunneling. Next, a gate pulse 
can then be applied such that the triplet-like level fi c ,T* is 
above the Fermi level of the left lead (c). A singlet-triplet 
spin-flip will lead to a tunneling out of an electron via this 
level. Tunneling events can be detected by charge sensing 
via an adjacent quantum point contact [7 j. These three 
pulse stages, a to c, are shown schematically in Fig. [5] 



Transport through the Singlet-like channel 

The strong line on the left edge of the singlet-like trian- 
gle in Fig. 1 of the main text can be explained by Fig. [6] 
At C, the fi c ,s* level of the left dot is resonant with the 
Fermi level of the left lead, /i^, as shown in Fig. [6^a). 
Also, /i£(2,o) an d Ms(i,i) are aligned with each other as 
shown in Fig. kmc). These two resonance conditions lead 
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(a) Step 1 



/'. .7 




l^d,T* 




Vd,S* 


1/^(2,0) 






^T(l.l) 


1^5(2,0) 


MS(1,1) 



Charging of left dot: system 
goes from ( I , I ) to S*(2, 1 ) state 



(b) Step 2 



A t T(2,0) 



MS(2,0) 



Hd,T* _ 



Mr(i,i) 



MS(1,1) 



Discharging the right dot: system 
goes from S*(2,l)to S(2,0) 



(c) Step 3 



Mr(2,o) 



MS(1,1) 



Interdot tunneling: system 
goes from S(2,0) to S( I , I ) 



Singlet-like channel transport at C 

FIG. 6: Diagrams of singlet-like channel transport at point 
C. Purple dots at the bottom of each diagram indicate the 
electron number at the beginning of each step in the transport 
cycle, (a) Starting from a (1,1) state, charging of the left dot 
from the /jl c ,s* level brings the system to the 5* (2, 1) state, 
(b) From the 5* (2, 1) state, the left dot ends up in the 5(2, 0) 
state as the right dot discharges from the /J,d,s* level, (c) 
Inter-dot tunneling from (2,0) to (1,1) completes the cycle. 



to a sharp current peak at C. Going away from C along 
the left edge of the triangle (Fig. 1, main text), the right 
dot falls in energy but the left dot maintains its reso- 
nance with the left lead, giving rise to a bright line along 
the left edge. In the direction where both dot energies 
fall in tandem, loading of the left dot into /jL c ,s* from 
the left lead is suppressed due to the exponential energy 
dependence in tunneling [T3] . 



Spin-flip Cotunneling in Transport Data at Forward 

Bias 



Here we discuss the interpretation of transport data 
taken at forward bias of Vsb — +0.526mV, originally re- 
ported in Ref. [13| and shown in Fig. [5] and show that 
certain features of the data arise because of spin-flip co- 
tunneling. 

In the three-electron regime, transport at forward bias 
cycles from (1,1) — » (2,0) — » (2,1). The bias triangles 
for forward bias are shown in Fig. [7^b). In this regime, 
transport is spin blockaded within the region bounded 
by jScruj of Fig. [7^b), due to the fact that /iT(i,i) is be- 
low the //t(2,o) l eve l as shown in Fig. [8]^c) . Any mecha- 
nism that allows spin relaxation to the 5(1, 1) state lifts 
the spin blockade, and typical relaxation mechanisms are 
expected to give rise to uniform current over the entire 
blockade region ^ftcruj. However, we observe a line of high 
current, shown by the dotted line in Fig. [7| which indi- 
cates that there is a preferentially high relaxation rate 
along this feature. We explain this line of strong current 
by a spin- flip cotunneling process. 

Point E (Fig. [7]) lies along the line of high current in the 
spin blockade region. At this point, the lower (solid line) 
of the two Hd,T* levels of the left dot are lined up with 
the right lead. This allows cotunneling of an electron 



from the right lead into the left dot to form the triplet- 
like T*(2,l) state. (Fig. |Jd)). Discharging of the left 
dot from the higher (dashed line) of the /jl^t* chemical 
potential results in the singlet 5(1, 1) state, which is not 
blockaded. The double dot is then able to complete the 
transport cycle until the next time the double dot gets 
spin blockaded. 

It is because the dotted line in Fig. [7] lies in the right 
dot energy axis that the spin blockade is lifted at a much 
faster rate. This is due to the fact that the left dot energy 
is aligned with the Fermi level of the right lead, allowing 
cotunneling to occur at a fast rate. As a result, spin 
relaxation happens at a much higher rate along this line. 
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FIG. 7: Transport data and schematic of bias triangles at forward bias, (a) Transport current /sd as a function of two gate 
voltages, Vg (V) and Vcs (V), taken at a forward bias source-drain voltage of Vsd = +0.526mV. These data were reported 
originally in Ref. [13] The feature shown by the dotted line is explained by the spin-flip cotunneling process, (b) Schematic 
of the bias triangles of transport data taken at the forward bias. Bias triangles corresponding to transport through ground 
(excited) states are shown in the red (blue) on the left (right). In the three electron regime (upper triangles), the dotted line 
lies in the region jSctcj, where the triplet T(l, 1) state is spin blockaded. Spin relaxation to the 5(1, 1) state lifts the spin 
blockade and allows transport to resume. This is possible within the region bounded by jSau and should give rise to a uniform 
current within the region. However, the presence of the line of high current (dotted) within 'jSauj suggests a preferentially high 
spin relaxation rate which is explained by the spin- flip cotunneling process shown below in Fig. [8] 



(a) 



Hd,s* 

MT(2,0) 



/AS(2,0) 



M 



Mr(i,i) 



VS(1,1) 



(b) 



'Hd,S* 

MT(2,0) 



^5(2,0) 



M 



Hc,S* 



f^T (1,1) 



VS(1,1) 



(c) 



Vd,S 
^T(2 ,0\ 
^S(2,0) 



(d) 




(e) 



1ms 

Mr(2,0) 
^5(2,0) 

99 



Vc,S* 



FIG. 8: Processes involved in the lifting of the triplet T(l, 1) spin blockade through cotunneling induced spin flip at point E. In 
the three-electron regime (upper triangles) of Fig.[7jb), transport is spin blockaded within the region bounded by jdauj. Purple 
dots at the bottom of each diagram indicate the electron number at the beginning of each step in the transport cycle, (a) 
Starting from any one of the two (2,0) states, the 5* (2, 1) or T*(2, 1) state is loaded from the right lead onto the left dot at the 
l^c,s* or /i c ,T* chemical potential respectively, (b) The left dot may discharge from the lower (solid line) of the [id, s* or /Jd,r* 
chemical potential, resulting in a T(l, 1) triplet. (Discharge from any one of the two dashed levels is also possible, but it results 
in a singlet 5(1, 1) that is not spin blockaded.) (c) The T(l, 1) triplet state is spin blockaded as the /Xt(i,i) chemical potential 
lies below the /xt(2,o) level, (d) Along the dashed line of Fig. u\ the lower (solid line) of the two /j,d,T* chemical potentials is 
lined up with the Fermi level of the right lead, allowing the loading of the excited (2,1) state through cotunneling from the right 
lead, (e) Subsequent discharge of a left dot electron could be from either one of the /id,T* levels, with discharge from the upper 
(dashed line) of the /J>d,r* levels resulting in the singlet 5(1, 1) state which can then complete the transport cycle by going into 
the 5(2,0) state through inter-dot tunneling. Transport thus resumes until the next triplet (1,1) spin blockade. The key is 
that from either the singlet-like or triplet-like (2,1) states, discharge from the left dot at the chemical potentials indicated by 
the dashed or solid lines results in a singlet or triplet (1,1) state respectively. The former allows transport to resume, whereas 
the latter results in a spin blockade. The dotted line of Fig. [7] lies in the direction where the left dot energy is fixed, leading to 
the line of high current. 



